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Improved modeling of Coulomb effects in 
nanoscale Schottky-barrier FETs 

Klaus Michael Indlekofer, Joachim Knoch, and Joerg Appenzeller 



Abstract — We employ a novel multi-conflgurational self- 
consistent Green's function approach (MCSCG) for the simula- 
tion of nanoscale Schottky-barrier field-effect transistors. This 
approach allows to calculate the electronic transport with a 
seamless transition from the single-electron regime to room 
temperature field-effect transistor operation. The particular im- 
provement of the MCSCG stems from a division of the channel 
system into a small subsystem of resonantly trapped states for 
which a many-body Fock space becomes feasible and a strongly 
coupled rest which can be treated adequately on a conventional 
mean-field level. The Fock space description allows for the 
calculation of few-electron Coulomb charging effects beyond 
mean-field. 

We compare a conventional Hartree non-equilibrium Green's 
function calculation with the results of the MCSCG approach. 
Using the MCSCG method Coulomb blockade effects are demon- 
strated at low temperatures while under strong nonequilibrium 
and room temperature conditions the Hartree approximation is 
retained. 
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I. Introduction 

ONE of the major challenges for the simulation of 
nanoscale field-effect transistors (FET) consists in an 
adequate description of the Coulomb interaction within the 
transistor channel: a proper simulation approach has to account 
for the Coulomb interaction of a few fluctuating electrons and 
at the same time has to be able to describe non-equilibrium 
transport in an open nanosystem. 

For the correct many-body description of the Coulomb 
interaction with the inclusion of contact coupling and nonequi- 
librium injection conditions, Fock space approaches such as 
realtime renormalization group [1], [2] or the Fock space 
Green's function [3] are available. This class of methods pro- 
vides kinetic equations in Fock space, taking renormalization 
and dissipation due to the coupling to the contacts (at least to 
some extent) and the Coulomb interaction into account. Since 
these approaches involve the 2 Ar -dimensional Fock space for 
the considered single-particle basis of N states, they typically 
are restricted to N < 10 for practical reasons. In the limit 
of small coupling, the Fock space description of the system 
can be approximated by a reduction to rate equations [4], 
[5], which deal with many -body eigenstates of the uncoupled 
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Hamiltonian obtained via exact diagonalization [6], [7], [8], 
[9], containing a full description of the Coulomb interaction. 

Realistic modeling of a ID semiconductor nanotransistor 
typically involves a number TV of single-particle states (orbital 
or sites) of up to a few hundred, rendering a full numeri- 
cal Fock space description impossible due to the exponen- 
tial scaling of the resulting many-body space dimensions. 
Furthermore, most of the potentially current-carrying single- 
particle states are stronlgy coupled to the contacts, and thus, 
the picture of a weakly coupled system in general becomes 
inadequate. The nonequilibrium Green's function (NEGF) 
approach [10], [11], [17], [18] in a mean-field approximation 
provides reasonble scalability, however, in principle lacks the 
description of few-electron Coulomb charging effects which 
become apparent in the case of resonantly trapped states in 
particular at lower temperatures. A possible solution is the 
combination of the numerically well scaling mean-field NEGF 
with a Fock space description for those states where many- 
body Coulomb effects may become important for the device 
characteristics. In this context, we have recently proposed a 
multi-configurational selfconsistent Green's function approach 
(MCSCG) [9], [12] for the realistic simulation of nanodevice 
systems under application-relevant conditions with reasonable 
numerical efforts. 

In the following sections, we will outline the main ideas 
and the algorithm behind the MCSCG and demonstrate its 
strengths by comparing a conventional Hartree NEGF calcu- 
lation with the results of the MCSCG approach, providing 
significantly extended information to our recent conference 
contribution [13]. Using the MCSCG, Coulomb blockade 
effects are demonstrated at low temperatures while under 
strong nonequilibrium and room temperature conditions the 
Hartree approximation is retained. 

II. MULTI-CONFIGURATIONAL APPROACH 

In order to handle systems with a large number N of single- 
particle states, the main idea of the MCSCG is to divide 
the system into two subsystems: Resonantly trapped (i.e. 
weakly coupled) states and those states that couple strongly 
to the contacts as depicted in Fig. Q Within the subspace of 
7Y' < 10 resonantly trapped states, which require a many-body 
description of the Coulomb interaction, a Fock space method 
will be applied, whereas the rest (N—N 1 ) is treated adequately 
on an approximated NEGF level. The eigenstates of 
the many-body statistical operator (or Hamiltonian, depending 
on the Fock space method) within the resonantly trapped 
subspace will be referred to as configurations with weights 



IEEE TRANSACTIONS ON ELECTRON DEVICES, VOL. I, NO. 1, NOVEMBER 2006 



2 



single particle space 
e.g. 50 sites 



Basis: 

Eigenstates of single 
particle density matrix 
"natural orbitals 



Fock space: 2 50 dim 




resonantly trapped: 2 6 = 64 



Fig. 1, Example for the division of the single-particle space into two 
subspaces. Only for resonantly trapped states (6 in the example), a Fock 
space description is employed. The rest (evanescent) is treated by a mean- 
field description. 



{w K } corresponding to the respective eigenvalues. Thus, the 
configurations and their weights follow from a Fock space 
calculation, taking the detailed Coulomb interaction within this 
subspace into account. The many-body statistical operator p 
of the considered system can be written in the general form 



st K 



(1) 



where P K = \k)(k\ denotes the projection operator cor- 
responding to the eigenstate and p re st is the many- 
body statistical operator of the rest, which may depend on 
the configuration k. Motivated by this form, we define a 
configuration-averaged Green's function G of the system as 



(2) 



where G[k) corresponds to P K rg> Prestin] and shall ful- 
fill Dyson's equation with a suitable contact coupling plus 
Coulomb selfenergy approximation E[«], In the simplest case, 
S[k] may be of a decoupled mean-field form [9], which 
is adequate for temperatures well above the corresponding 
Kondo temperature. 

As for the Fock subspace of the resonantly trapped states, 
the projected many-body Hamiltonian contains the reduced 
single-particle and Coulomb terms. Furthermore, coupling to 
the contacts with nonequilibrium carrier injection and cou- 
pling (tunneling) to the rest of the system is described by 
means of selfenergy kernels, depending on the chosen Fock 
space method. For the latter, various choices are possible, 
for example: Exact diagonaliztion with Dyson's equation as 
subsidiary condition [9], [12], real-time renormalization group 
(RTRG) [1], [2] or Fock space Green's functions [3]. In 
the following, we will discuss the first option, based on 
exact diagonalization. In this case, the many-body statistical 
operator is assumed to be diagonal in the eigenbasis of the 
resonant subspace Hamiltonian. Here, the eigenvalues w K are 
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Fig. 2. Flowchart of the MCSCG algorithm. In the shown case, exact 
diagonalization is employed as the simplest Fock subspace method. Arrows 
visualize the flow of the selfconsistency loop. 



determined such that the resulting many-body Green's function 
Gmb in Lehmann representation fulfills Dyson's equation 
within the resonant subspace. In the simplest implementation, 
{w K } is chosen such that the spectral peaks of G MB match 
those of G< [9]. 

As an overall selfconsistency condition, the resonantly 
trapped states experience a mean-field interaction of the rest, 
whereas the rest is subject to the set of selfenergies 

originating from the resonant many-body configurations n and 
its own mean-field interaction. Finally, for the identification of 
resonantly trapped states, the single-particle eigenstates of the 
single-particle density-matrix are employed (so-called natural 
orbitals). The latter follows directly from the Green's function 
G < as part of the multi-configurational selfconsistency proce- 
dure. (Note that each individual G[n] need not be selfconsis- 
tent with its respective S[k].) In turn, resonantly trapped states 
are defined as those single-particle eigenstates that exhibit a 
level broadening (determined from the selfenergies) below a 
given threshold. Fig. [2] illustrates the details of the algorithm 
as a flowchart. 

For the calculation of expectation values of single-particle 
observables (e.g., electron density, current, spin density, etc.), 
the selfconsistent Green's function G K is employed as an 
approximation for the unknown exact G < , Optionally, within 
the resonant subspace, one can use the many-body result (such 
as Gmb, the reduced statistical operator, dissipation kernel, 
etc.) for the evaluation of arbitrary expectation values, in 
particular contour-ordered correlation functions of arbitrary 
order. 

III. Simulation of Nanoscale Schottky-barrier 
FETs 

In order to demonstrate the strengths of the MCSCG ap- 
proach, we consider a one-dimensional (ID) coaxially gated 
nanowire transistor with Schottky-barrier source and drain 
contacts since deviations from a mean-field approximation 
become most apparent in a system with quasi-bound states. 
Fig.|3]shows a schematics of such a nanowire transistor where 
we assume a channel length of L = 20nm, a diameter of 
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Fig. 3. Schematic sketch of a nanowire MOSFET with coaxial gate. 



d c h = 4nm and a gate oxide thickness of d ox = lOnm. (Since 
the spin degree of freedom is included and a site spacing of 
a = lnm is employed, we thus consider N = 40 single- 
particle states.) Such a nanowire transistor can in principle be 
realized with a semiconducting InGaAs nanowire with SW2 
as gate dielectric. It has been shown that the electrostatics of 
coaxially gated nanowire transistors can be well described by a 
modified one-dimensional Poisson equation [16]. This Poisson 
equation allows to easily calculate the Coulomb Green's 
function which in turn enables the description of the classical 
electrostatics and the screened interaction between electrons 
on equal footing [9]. For the following simulation results, 
the simplest MCSCG variant has been employed, based on 
exact diagonalization (with N' — 6 resonantly trapped states 
yielding 64 Fock space dimensions) with a decoupled static 
selfenergy form. 

Fig. |4] visualizes the simulated drain current Id for the 
single-electron transport regime (T = 77K) as a color plot. 
In contrast to the Hartree-only calculation (Fig. Ela)), the 
MCSCG approach (Fig. EJb)) correctly reveals diamond-like 
shaped patterns due to the quantized Coulomb interaction (as 
predicted by the orthodox theory and observed in experi- 
ments). While the MCSCG treatment is able to cope with the 
mixture of many-body configurations, the Hartree theory only 
provides a mean interaction potential for the description of the 
Coulomb interaction. 

In addition, Fig. |5] shows Id(Vgs) curves for different 
drain voltages Vbs- in me MCSCG case (Fig. |5] (b)), single- 
electron transport can be identified in terms of Coulomb 
oscillations for the two lowest Vbs> whereas the Hartree-only 
simulation (Fig. 0a)) lacks these features; the Hartree-only 
case exhibits broader peaks solely due to the single-particle 
levels of the system. However, with increasing Vds> both 
approaches become equivalent. (Note that the sub-threshold 
regime shows the regular behavior and has been omitted here.) 

Finally, Fig. [6] shows the room temperature (T = 300K) 
characteristics. Apart from the slight modulation in the MC- 
SCG calculation (Fig. |6jb)), which is a remnant of the 
Coulomb oscillation, the Hartree (Fig. |6ja)) and MCSCG 
(Fig. (b)) results are in good agreement. Effects beyond 
a mean-field picture of the system will obviously have a 
significant impact on application-relevant device properties 
such as the system capacitance and the transconductance. 
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Fig. 4. (a) Color plot of the total drain current for the Hartree-only case at 
T = 77K (blue = OnA, red = 5nA). (b) Color plot of the total drain current 
for the MCSCG case at T = 77K. (blue = OnA, red = 5nA). 



IV. Conclusion 

In summary, we have compared for the first time the conven- 
tional Hartree NEGF with the MCSCG and have shown that 
the multi-configurational approach is able to describe single- 
electron charging effects in the low temperature limit for a 
realistic FET structure. In case of strong nonequilibrium (with 
an almost depleted channel) and room temperature conditions, 
the MCSCG and the well-established Hartree approximation 
lead to equivalent results for the discussed example of a 
nanowire MOSFET. As such, the MCSCG yields a seamless 
transition from the single-electron transport regime to transis- 
tor operation at room temperature. For realistic FETs with a 
large number of sites where a full Fock space formulation 
becomes impossible, the MCSCG permits a selfconsistent 
Fock space treatment of states which are responsible for few- 
electron charging effects. 
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